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O ' Abstract 

H ' 

! We present the properties of Lya emitters (LAEs) at z = 5.7 in the Subaru Deep Field. A photometric 

sample of 89 LAE candidates is constructed from narrow-band (NB816) data down to A5816 = 26.0 (AB) 
in a continuous 725 arcmin 2 area. Spectra of 39 objects satisfying the photometric selection criteria for 
LAEs were obtained with Subaru and Keck II Telescopes, among which 28 were confirmed LAEs, one 
was a nearby galaxy, and eight were unclassified. We also obtained spectra of another 24 NB816-cxcess 
J~j \ objects in the field, identifying six additional LAEs. We find that the Lya luminosity function derived 



from the photometric sample is reproduced well by a Schechter function with L* = 7.9^2 x 10 42 er § s_1 
and <f>* = 6.3^2 q x 10 -4 Mpc~ 3 for a = —1.5 (fixed) over the whole luminosity range of L ~ 3 x 10 42 - 
3 x 10 43 erg s _1 . We then measure rest-frame Lya equivalent widths (EWs) for the confirmed LAEs, to 
find that the median among the 28 objects satisfying the photometric selection criteria is Wq = 233 A. We 
infer that 30% - 40% of LAEs at z = 5.7 exceed W [ = 240 A. These large-EW objects probably cannot 
be accounted for by ordinary star- forming populations with a Salpeter IMF. We also find that LAEs with 
fainter far-UV luminosities have larger EWs. Finally, we derive the far-UV luminosity function of LAEs 
down to Mtjv — —19.6 using the photometric sample, and compare it with that of Lyman-brcak galaxies 
(LBGs). We find that as high as about 80% of LBGs at z ~ 6 have Wq > 100 A, in sharp contrast to 
lower-z counterparts. 

Key words: galaxies: evolution — galaxies: high-redshift — galaxies: luminosity function, mass 
function — galaxies: photometry 
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1. Introduction 

Lya emitters (LAEs) are commonly seen in the high- 
redshift universe. Most of them are thought to be very 
young star-forming galaxies, although some may have an 
active nucleus and some could be primeval galaxies in a 
gas cooling phase before onset of the initial star formation 
or in superwind phases just after the initial star formation. 
The high number density of LAEs, roughly comparable 
to that of Lyman- break galaxies (LBGs), suggests that 
they arc important in understanding the early evolution 
of galaxies. 

LAEs at z > 5 arc of particular interest in galaxy for- 
mation and cosmology, since they are not only expected 
to be extremely young but also an invaluable probe of the 
universe during or just after reionization. They can be 
used to infer the reionization epoch of the universe via 
the change in their number density (Haiman & Spaans 
1999; Rhoads & Malhotra 2001; Haiman 2002; Malhotra 
& Rhoads 2004; Stern et al. 2005; Kashikawaet al. 2005). 
However, although there has been rapid progress in the 
observation of LAEs at z > 5 (e.g., Rhoads & Malhotra 
2001; Ajiki et al. 2003, 2004; Kodaira et al. 2003; Hu et 
al. 2002, 2004; Maier et al. 2003; Rhoads et al. 2003; 
Ouchi ct al. 2005; Taniguchi et al. 2005), the data ob- 
tained to date are not sufficient to place strong constraints 
on their fundamental properties such as the luminosity 
function (LF), in contrast to the data of LBGs at simi- 
lar rcdshifts. Usually, narrow-band surveys are used to 
detect LAEs. Their apparent faintness in both narrow 
bands (i.e., Lya) and broad bands (off Lya), as well as a 
small survey volume defined by the wavelength coverage 
of the narrow-band filter times the survey area, makes it 
difficult to obtain a large sample. Samples from several 
small surveys by different workers could be combined to 
form a larger sample. However, since the selection criteria 
of LAEs and the survey depth differ among the surveys, 
LFs calculated from such combined samples may have sys- 
tematic uncertainties difficult to control. 

In this paper, we report on the results of a very deep 
survey of LAEs at z = 5.7 in a wide, continuous area of 725 
arcmin 2 . The survey uses deep imaging data in six optical 
bandpasses including the narrow band NB816 (A c = 8150 
A and FWHM = 120 A) targeted at the Lya emission at 
z ~ 5.7 ± 0.05. All of the imaging data were taken with 
the Suprime-Cam (Miyazaki et al. 2002) at the Subaru 
Telescope (lye et al. 2004). From these data, we construct 
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a sample of 89 photometrically selected LAEs down to 
iVB816 = 26.0. By followup spectroscopy with the Subaru 
and Keck II Telescopes, we obtained a sample of 34 LAEs 
with confirmed rcdshifts. 

Combining these samples with Monte Carlo simula- 
tions, we calculate the Lya LF, which is the most reli- 
able measurement for the LF of z = 5.7 LAEs published 
to date. This LF serves as a 'zero point' for the LF in 
the fully ionized universe (z < 6), which is to be compared 
with the LF at z > 6 to probe the reionization history of 
the universe. A companion paper (Kashikawa et al. 2005) 
compares the LF at z = 5.7 obtained in this paper with 
that for z = 6.5 LAEs, and discusses the ionization state 
of the intcrgalactic medium at z = 6.5. 

We also derive accurately the far UV continuum LF and 
the distribution of Lya equivalent widths (EWs), both of 
which have not yet been strongly constrained for LAEs at 
z > 5. Comparing the far UV LF of our LAEs with that 
of LBGs at z ~ 6 given in the literature, we examine what 
fraction of LBGs have EWs large enough to be detected 
as LAEs, and discuss evolution of Lya properties of star- 
forming galaxies. 

The plan of this paper is as follows. Section 2 describes 
the photometric and spectroscopic data used in this study. 
The confirmed sample of 34 LAEs is constructed in §3. 
The photometric selection of LAEs is made in §4. Section 
5 is devoted to results and discussion, and a summary 
is given in §6. AB magnitudes are used in this paper. 
We assume a flat universe with = 0.3, J7a = 0.7, and 
H = 70 km s^ 1 Mpc" 1 . 

2. Data 

2.1. Imaging Data 

The SDF has deep, public Subaru/Suprime-Cam data 
in seven bandpasses, B, V, R, i' , z', NB816, and NB921, 
obtained for the Subaru Deep Field Project (Kashikawa 
et al. 2004; see also Maihara et al. 2001 for the Subaru 
Deep Field). This project is a large program of Subaru 
Observatory to carry out a deep galaxy survey in a blank 
field in optical and NIR wavelengths to study distant 
galaxies. B, V, R, i' , and z' arc standard Johnson and 
SDSS broad bands, and NB816 and NB921 are narrow 
bands whose central wavelength and FWHM are (8150 A, 
120 A) and (9196 A, 132 A), respectively. In this paper, 
we use the B, V, R, i', z', and NB816 data to study LAEs 
at z = 5.7. The exposure times and the 3er limiting mag- 
nitudes on a 2" aperture are: 595 min and 28.45 mag (B), 
340 min and 27.74 mag (V), 600 min and 27.80 mag (R), 
801 min and 27.43 mag («'), 504 min and 26.62 mag \z'), 
600 min and 26.63 mag (NB816). These images have been 
convolved to a seeing of 0."98 (FWHM) and an identical 
sky coverage of 29/7 x 36/7. The pixel scale of the images 
is 0."202 pixel -1 . The selection of z = 5.7 LAEs is made 
in the R — z' vs i' — NB816 plane, and the B and V data 
are also used to further remove foreground objects. 

Object detection and photometry were made using 
SExtractor version 2.1.6 (Bcrtin & Arnouts 1996) on all 
six images. The NB816-band image was chosen to de- 
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tect objects. If more than 5 pixels whose counts were 
above the 2 a s ky were connected, they were regarded as 
an object. In total, 82,212 objects were detected down to 
the 3<T limiting magnitude of NB816. For each object de- 
tected in the NB816 image, a 2"-diameter magnitude was 
measured for each passband to derive the colors of each 
object. We adopted MAG_AUTO for the total NB816 
magnitude. All the observed magnitudes and colors are 
corrected for Galactic absorption using the dust maps of 
Schlegel, Finkbeiner, & Davis (1998). 

2.2. Spectra with FOCAS and DEIMOS 

We made spectroscopic observations of 63 objects se- 
lected from the NB816-detected catalog with FOCAS 
(Kashikawa et al. 2002) on the Subaru Telescope and 
DEIMOS (Faber et al. 2003) on the Keck II Telescope. All 
63 objects are NB816-excess objects with i' — NB816 > 1.0 
and NB816 < 26.0. We first assigned slitlets to objects ei- 
ther with R — z' > 1.0 or undetected in R; they are likely 
LAE candidates. We then included objects irrespective 
of R — z' color in the target list if slitlets were available, 
in order to evaluate the completeness and contamination 
of the color-selected sample of LAE candidates to be con- 
structed in Section 4. 

2.2.1. FOCAS Spectroscopy 

A total of 39 objects were observed with FOCAS in 
2002 and 2004 in the multi-slit mode. Among them, three 
objects were observed on UT 2002 June 7 and 9 with a 
300 lines mm -1 grating and an 058 order-cut filter using 
two masks. The typical spectral coverage was 5400 A - 
10000 A with a pixel scale of 0.10 A pixer 1 . The use of 
0."8 slits gave a resolution of 9.5 A at 8150 A. The spatial 
resolution was 0."3 pixel -1 by 3-pixel on-chip binning. 
The integration time per mask was 12600 seconds. The 
sky condition was good with a seeing size of 0."4 - l."0. 
We also obtained spectra of the standard stars Hz 44 and 
Feige 34 for flux calibration. The data were reduced in a 
standard manner. 

The remaining 36 objects were observed on UT 2004 
April 24-27 with the same type of grating and order-cut 
filter. Six masks were used in total. The use of 0."6 
slits gave a resolution of 7.1 A at 8150 A. The spatial 
resolution was 0."3 pixel" 1 by 3-pixel on-chip binning. 
The integration time per mask was 12000-16800 seconds. 
The sky condition was good with a seeing size of 0."4 - 

0. "8. The standard stars Hz 44 and Feige 34 were observed 
for flux calibration. The data were reduced in a standard 
manner. 

2.2.2. DEIMOS Spectroscopy 

Twenty-four objects among the 63 were observed with 
DEIMOS on UT 2004 April 23 and 24. We used four 
MOS masks with an 830 line mm" 1 grating and a GG495 
order-cut filter. The spectral coverage was - 5000 A - 
10000 A, with a central wavelength of 7500 A for one MOS 
mask and 8100 A for the other three. The slit width was 

1. "0, giving a spectral resolution of 3.97 A. The integra- 
tion times of individual masks were 7000 - 9000 seconds, 



and the typical seeing sizes were 0."55 - l."0. Spectra 
of the standard stars BD+28d4211 and Feige 110 were 
taken for flux calibration. The data were reduced with 
the spec2d pipeline. 1 

3. Spectroscopic Sample of z — 5.7 LAEs 

3.1. Classification of Spectra by Weighted Skewness 

We classify the 63 objects with spectroscopic observa- 
tions into LAEs and nearby emission-line objects using 
the features of emission lines. NB816-excess objects can 
be either nearby emitters (Ha, or [Oiii], or [On] emitters), 
or LAEs at z ~ 5.7. Typically, the Lya emission lines of 
high-rcdshift galaxies are single, asymmetric lines, with 
the profile blueward of 1216 A erased almost completely 
by neutral hydrogen gas. Thus, they can be discriminated 
from nearby emission lines, which are either single, sym- 
metric lines, or doublet lines. 

Objects with doublet emission lines are easy to clas- 
sify; they are either [On] or [Oiii] emitters depending on 
the separation of the doublet lines. Objects with single 
emission lines are either LAEs, or Ha emitters, or [On] 
emitters with unresolved doublet lines. To classify sin- 
gle emission-line objects, we introduce weighted skewness, 
Sw, which quantifies the asymmetry of emission lines. A 
similar quantity (skewness) has been applied by Kurk et 
al. (2004) to identify an LAE at z = 6.518. 

For a given one-dimensional spectrum of an emission 
line, f(xi) (i = l,2,...,n), where Xi is the i-th pixel along 
wavelength and f(xi) is the flux in the z-th pixel, the 
weighted skewness is defined as 

1 - 

S W = (AlO.r - Aio.b) x j-^ }Xxj - x) 3 f(xi). (1) 

i=l 

Here, I = J2 fi x i)'> i = 7E Xif(xi); & = 
\JlJ2( x t -x) 2 f( x i)'> Aio, r and \ w , b are the wave- 
lengths where the flux drops to 10% of its peak value 
at the red and blue sides of the emission, respectively. 
In the above formula, (pELif 1 ' - ^) 3 /! 1 ') measures 
the skewness of the line, and Xio, r — Aio,& measures the 
width of the line. Since the Lya emission of high-redshift 
galaxies tends to be wider than other emission lines 
of nearby galaxies in the observed frame, we include 
the factor Aio, r — Aio,& in Sw to enhance the difference 
between Lya and other lines. Sw is larger for objects 
with higher asymmetries and/or larger line widths. Note 
that for an identical object, spectroscopy with a higher 
dispersion generally gives a larger Sw value. A full 
description of Sw is found in Kashikawa et al. (2005). 

Lya emission lines at high redshifts typically have large 
positive Sw values. On the other hand, [Oiii] and Ha 
lines are nearly symmetric, i.e., Sw — 0. The Sw of [On] 
emitters is also expected to be small (Sw is measured for 
the doublet as a whole); for all resolved [On] emitters in 
our sample, the A3726 line is weaker than the A3729 line, 

This pipeline was developed at UC Berkeley with support from 
NSF grant AST-0071048. 
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Fig. 1. Weighted skewness of the 63 objects with spectro- 
scopic observations. The red circles, blue triangles, and green 
crosses indicate LAEs, nearby objects, and unclear objects, 
respectively. The objects with open blue circles have multiple 
emission lines. Ordinate shows wavelength ratio, a\. 

meaning a negative Sw (See also Rhoads et al. 2003). 
This means that Sw can discriminate LAEs from [On] 
emitters even when the [On] doublet is not resolved in a 
spectrum due to a low spectral resolution and/or a low 
signal-to-noise ratio. 

Figure 1 plots Sw for the 63 objects. The ordinate 
represents, for a reference, the wavelength ratio (a\) de- 
fined by Rhoads et al. (2003); a\ is defined as a\ = 
(Aio,r — Ap)/(A p — Aio.fi), where X p is the wavelength of 
the peak flux density. A loose correlation is found be- 
tween a\ and Sw hi Fig- 1- Objects marked by a large 
open circle are nearby emitters with definite classification 
based on their emission lines; they are either Ha emitters 
(from the presence of the [Oin] line or from a very flat 
and strong continuum over the whole wavelength range), 
or [Oni] emitters (from doublet lines at 4959 A and 5007 
A), or [On] emitters (from doublet lines at 3726 A and 
3729 A). All these nearby objects have Sw < 3. A visual 
inspection of the 63 spectra reveals that most of the ob- 
jects with Sw > 3 have an asymmetric emission-line shape 
typical of high-redshift Lya emission. 

On the basis of these results, we classify the spectro- 
scopic objects into three classes as follows. We primarily 
use Sw, but in some cases we also use B and V photom- 



etry as additional information, since LAEs at z = 5.7 are 
essentially invisible in our B and V data due to strong 
absorption by the IGM. If an object has the definite fea- 
ture of the [On] or [Oin] or Ha line, it is a nearby object. 
If an object has a single emission-line with Sw > 3, it is 
regarded as an LAE (irrespective of the B and V fluxes); 
we assume that those with detected B or V fluxes are 
contaminated by a faint foreground object. If, on the 
other hand, a single-line object with Sw < 3 is detected 
in either B or V (i.e., brighter than the 2a limit in the 
corresponding bandpass), then it is classified as a nearby 
object. Objects which do not satisfy either of the above 
three categories are classified as 'unclear' objects; unclear 
objects are single-line emitters with Sw < 3 but invisible 
in both B and V. Among the 63 objects, 34, 19, and 10 
are classified as LAEs, nearby objects, and unclear ob- 
jects. Among the 34 LAEs, six are detected in either B or 
V. This rather high contamination rate, 6/34 = 18%, ap- 
pears to be accounted for by the observed surface number 
density of faint B/V sources. We find that the fraction of 
the sky occupied by faint (m > 26) B/V sources is about 
10%. Furthermore, the probability of a faint B/V source 
being within 1" of a random position in the sky is found to 
be about 20%. The contamination rate obtained for the 
spectroscopic sample may be compared to these rough es- 
timates of chance projection. 

Table 1 summarizes the photometry and spectroscopy 
of the 63 objects. Figures 2-4 show the spectra of the 
LAEs, nearby objects, and unclear objects, respectively. 
The luminosity functions and star formation properties of 
the line-emitting galaxies at z < 1.2 identified in our spec- 
troscopy will be discussed in Ly et al. (2006, in prepara- 
tion) . 

We infer that some portion of the unclear objects may 
be LAEs in reality This is because (i) not all LAEs may 
have as large asymmetry as the confirmed LAEs and (ii) 
since the resolution of the FOCAS spectra is significantly 
poorer than that of the DEIMOS spectra, the absolute 
value of Sw of an object based on a FOCAS spectrum will 
be smaller than that based on a DEIMOS spectrum. In 
this sense, the criterion of Sw > 3 is conservative in terms 
of the selection of LAEs. Among the ten unclear objects, 
eight are FOCAS objects. When estimating the selection 
completeness of the photometric sample of z = 5.7 LAEs 
in Section 4, we take account of the extreme possibility 
that all of the unclear objects are actually LAEs. 

Figure 5 shows the distributions of emission-line wave- 
lengths (the central value if the emission is a doublet) for 
all of our objects. The thick solid histogram in panel (a) 
represents the distribution for the 34 LAEs, and panel (b) 
plots the distributions for nearby objects (thin solid his- 
togram) and unclear objects (dotted histogram). In both 
panels, the dashed curve corresponds to the response func- 
tion of the NB816 band. It is found that while the dis- 
tributions of nearby and unclear objects are statistically 
consistent with being symmetric around the central wave- 
length of the bandpass, the LAEs' distribution is skewed 
toward shorter wavelengths. A natural explanation for 
this skewed distribution is as follows. LAEs with higher 
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Fig. 2. Spectra of LAEs. The ID number (in the 
NB816-detected catalog) and redshift are shown in the upper 
right corner of each panel. 'F' ('D') in parentheses indicates 
that the object was observed with FOCAS (DEIMOS). The 
vertical dotted lines indicate the center of the (single) emis- 
sion line. The scale on the y axis is marked in 0.5 X 10 — 18 erg 
s~ 1 cm~ 2 A -1 ; for panels in which a factor is shown in the 
upper left corner, multiply the scale by this factor to obtain 
a correct scale. 

redshifts are on average fainter in NB816 (even when the 
Lya luminosity is fixed) because of a lower fraction of the 
continuum emission rcdward of A res t = 1216 A falling into 
the NB816 bandpass. Hence they tend to be excluded 
from a list of spectroscopic targets. 

The dotted curve in panel (a) represents the distribu- 
tion of modeled LAEs with iVi?816 < 26.0 in a mock sam- 
ple; this sample is constructed based on the same selec- 
tion criteria as for the observed sample (see the next sec- 
tion) from a large catalog of LAEs generated by Monte 
Carlo simulations for a range of LAE equivalent widths 
and luminosity function parameters (see Section 5). In the 
simulations, LAEs are distributed uniformly in comoving 
space over the redshift range 5.6 < z < 5.8, which corre- 
sponds to an observed Lya wavelength of 8020 to 8270 A. 
The thick histogram is found to agree fairly well with the 
dotted curve, suggesting that prominent large-scale struc- 
tures along the line of sight are not present in our LAE 
sample. 

Hu et al. (2004) also found a skewed distribution in 
their z = 5.7 LAEs selected from NB816 data in the SSA22 
field, and argued that the true center in wavelength of the 
NB816 filter may be shifted toward shorter wavelengths 
by about 25 A than the measured one. Although we can- 



Fig. 2. Fig. 2 - continued. The second bottom and bot- 
tom panels in right-hand side show typical spectra of the 
sky background obtained with FOCAS and DEIMOS, respec- 
tively. The normalization is arbitrary. 

not completely rule out this possibility, the wavelength 
distribution of our sample is explained without assuming 
such a shift. 

3.2. Confirmed Sample of z = 5.7 LAEs 

We have 34 LAEs with spectroscopic confirmation. 
Their distribution in the i' — NB816 vs NB816 plane 
is slightly biased toward brighter NB816 magnitudes and 
redder i' — NB816 colors with respect to the photomet- 
ric sample, except for three objects of i' — NB816 ~ 1.2. 
Among the 34 objects, the faintest in NB816 is 25.89 mag. 

We measure the flux and FWHM of the Lya line from 
the spectrum of each object. In this paper, however, we 
adopt for Lya fluxes the values calculated from the imag- 
ing data by subtracting the continuum flux measured from 
the z'-band image from the flux falling into the NB816 
band, assuming /„ = constant; spectroscopic redshifts are 
used to determine the precise position of the Lya line in 
the NB816 response function. The reason for not using 
the Lya fluxes directly is that our multi-slit spectroscopy 
probably lost a non-negligible amount of Lya photons, 
because of the relatively narrow slitlcts (0."6 and 0."8 
for FOCAS and l."0 for DEIMOS) compared with the 
typical PSF sizes, and because the Lya emission could 
be extended spatially. Figure 6 shows Lya luminosities 
measured from the spectra against those from photom- 
etry (See Subsection 5.2 for the calculation of Lya lu- 
minosities based on the photometric data). While lit- 
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Fig. 3. Spectra of nearby objects. The ID number (in the 
NB816-detected catalog) and redshift are shown in the upper 
right corner of each panel. 'F' ('D') in parentheses indicates 
that the object was observed with FOCAS (DEIMOS). The 
vertical dotted lines indicate the identified emission lines (Ha 
or [Om] or [On]). Identification are unsuccessful for #28247, 
#29275, #34775, and #168136; these are classified as nearby 
objects because of a small skewness (Sw < 3) an d a bright 
B or V magnitude. The scale on the y axis is marked in 
1 X 10" 18 erg s -1 cm" 2 A -1 ; for panels in which a factor 
is shown in the upper left corner, multiply the scale by this 
factor to obtain a correct scale. 



tie systematic difference is seen between the luminosities 
from the DEIMOS spectra and from the photometry, the 
FOCAS spectra give luminosities only about one third 
those measured from the photometric data. Among the 
FOCAS objects, three were observed with a 0.8" slitlct, 
and their i(Lya(spec))/i(Lya(phot)) ratio is 70% on av- 
erage. These results are qualitatively consistent with nar- 
rower slitlets losing a larger fraction of Lya photons. It is 
difficult to correct L (Lya (spec)) for this effect, since the 
seeing size varied during the exposures and since different 
objects have different intrinsic sizes. For the same reason, 
we also measure the equivalent widths of the 34 objects 
from the photometric data, as described in Subsection 5.2. 

4. Photometric Sample of z — 5.7 LAEs 

We construct a photometrically selected LAE sample 
from the imaging data. Objects located in low-quality 
regions of the images are excluded from the selection pro- 
cess to obtain an LAE sample of uniform detection. The 
low-quality regions include areas around very bright stars, 
the four edges of the images, and the south-east region 
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Fig. 4. Spectra of unclear objects. The ID number (in the 
NB816-detected catalog) is shown in the upper right corner of 
each panel. 'F' ('D') in parentheses indicates that the object 
was observed with FOCAS (DEIMOS). The vertical dotted 
lines indicate the position of the single emission line. The 
scale on the y axis is marked in 0.5 X 10~ 18 erg s _1 cm -2 
A -1 ; for panels in which a factor is shown in the upper left 
corner, multiply the scale by this factor to obtain a correct 
scale. 



which corresponds to the field of view of the CCD whose 
quantum efficiency is only about two thirds those of the 
other nine CCDs. The effective area is thus reduced to 
725 arcmin 2 . 

Figure 7 plots i' — A^£?816 color against NB816 mag- 
nitude for all objects with ^5816 < 26.0 (N = 61,919). 
The vast majority of the objects are distributed around 
i' - NB816 ~ 0.2; their i' - NB816 distribution spreads as 
NB816 magnitude goes fainter, mainly due to the increase 
in photometric errors. Objects with large i' — ^5816 col- 
ors are emission-line objects (or some kind of late-type 
stars with i' — NB816 > 1), part of which arc LAEs at 
z = 5.7. The large symbols indicate objects with spec- 
troscopic observations: red circles, blue triangles, and 
green crosses indicate LAEs (N = 34), nearby objects 
(19), and unclear objects (10), respectively. For LAEs, 
open circles mean R magnitudes fainter than the 2a mag- 
nitude (N = 29). The magenta line shows the thresh- 
old of i' — NB816 adopted for the selection of LAEs: 
i' — NB816 > 1.5 (see below). The dotted line implies 
the i'-band 2a magnitude. The i'-band magnitudes of ob- 
jects undetected in i' (i.e., fainter than the 2a magnitude) 
have been replaced with the 2a magnitude. 

Figure 8 shows the distribution of objects in the R— z' 
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Fig. 5. Distributions of emission line wavelengths for objects 
with spectroscopic observations. Panel (a). The distribution 
for the 34 LAEs is shown by the thick solid histogram. Panel 
(b). The distributions for the 19 nearby objects (thin solid 
histogram) and the 10 unclear objects (dotted histogram). 
The dotted curve in panel (a) represents the distribution of 
model LAEs with 7VB816 < 26.0 in a mock sample (See text). 
In both panels, the dashed curve corresponds to the response 
function of the NB816 band. 

vs i' — NB816 plane. In this figure, if an object is un- 
detected (i.e., fainter than the 2a magnitude) in a given 
bandpass, its magnitude in that bandpass is replaced with 
the 2a magnitude. Objects on the R — z' = 1.16 line are 
those undetected in both R and z'. More than half (18/34) 
of the LAEs are invisible in z'. The red circles represent 
the spectroscopically confirmed LAEs (N = 34) ; open cir- 
cles show that they are undetected in R (N = 29). The 
blue triangles and green crosses indicate nearby objects 
(19) and unclear objects (10), respectively. The yellow 
stars correspond to the 175 Galactic stars given in Gunn 
& Stryker's (1983) spectrophotometric atlas. The spectro- 
scopically confirmed LAEs are found to be separate fairly 
well from the nearby objects. On the other hand, the un- 
clear objects are located in a similar area to the LAEs; 
this may suggest that a significant fraction of the unclear 
objects are LAEs in reality. 

The black solid curves show the tracks of model LAEs 
over a redshift range of 5.5 < z < 5.9 with five different 
intrinsic EW values of Lya in the rest frame: Wq = 0, 
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Fig. 6. Lya luminosities measured from the spectra are plot- 
ted against those based on the photometric data for the 34 
LAEs with spectroscopic confirmation. The filled and open 
circles indicate objects observed with DEIMOS and FOCAS, 
respectively. Three FOCAS objects observed through 0."8 
slitlcts arc marked with large open circles. 

20, 50, 150, and 300 A from left to right; filled squares 
correspond to z = 5.7. In this paper, we assume that the 
intrinsic profile of Lya emission is symmetric centered on 
zero velocity, and that foreground Hi gas reduces intrinsic 
EWs 50% (blue half). Thus, we set the rest-frame ap- 
parent equivalent width after truncation by Hi gas to be 
Wq = 0.5 x Wq] the quantity measured from the data is 
Wq. Note that the assumption of 50% reduction has not 
been verified for LAEs at z ~ 5.7 and may be an over- 
simplification. As stated in Dawson et al. (2004), the 
dust content and detailed kinematics of the galaxy ISM 
appear to play a significant role in determining the emer- 
gent Lya profile (e.g., Kunth et al. 1998; Stern & Spinrad 
1999; Mas-Hesse et al. 2003; Shapley et al. 2003; Ahn 
2004). LAEs at z ~ 5.7 may have a wide variety in the 
total absorption fraction of Lya photons. 

To calculate the model tracks, we construct the contin- 
uum spectral energy distribution of model galaxies using 
Kodama & Arimoto (1997) stellar population synthesis 
model, adopting a Salpeter IMF, a 10 7 yr age, and con- 
stant star formation. Note that LAEs without Lya emis- 
sion can have as red i' — NB816 colors as ~ 1.2. 

Based on this figure, we select z = 5.7 LAEs, by im- 
posing the following four criteria on the NB816-detected 
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NB816 

Fig. 7. i' - TVB816 plotted against TVB816 for objects with 
TV 73816 < 26.0. The large symbols indicate objects with spec- 
troscopic observations: red circles, blue triangles, and green 
crosses indicate LAEs (TV = 34), nearby objects (19), and un- 
clear objects (10), respectively. For LAEs, open circles mean 
R magnitudes fainter than the 2a magnitude. The magenta 
line shows the threshold of i' — TV 73816 adopted for the selec- 
tion of LAEs: i' - TV 73816 > 1.5 (see text). The dotted line 
implies the i'-band 2a limit. 

multi-color catalog: 

7V73816 < 26.0 (i) 

i' - NB816 > 1.5 (ii) 

{R — z' > 1.5 and mR < mn(2a)} or {ttir > ttir(2(t)} (iii) 

bib > mB(2cr) and my > mv(2cr) (iv) 

Objects satisfying these four criteria simultaneously are 
regarded as LAE candidates. The magenta solid lines in 
the figure indicate the selection boundary (in the case of 
tor < mn(2(j)) in the R — z' vs i' — NB81 6 plane. In 
total 89 LAE candidates are selected, which we call the 
photometric sample. 

Table 2 summarizes a breakdown of the spectroscopic 
sample (./V = 63) in terms of the selection criteria for our 
photometric sample. Thirty-nine out of the 63 objects 
satisfy the photometric selection criteria for LAEs; 28 are 
confirmed LAEs, one is a nearby galaxy, and ten arc un- 
classified. In the table, numbers in parentheses are the 
number of objects lying in the 'high-quality' region of 
725 arcmin 2 used to construct the photometric sample 
of LAEs. Among the 28 confirmed LAEs, one is out- 
side the 'high-quality' region, and thus is not included 
in the photometric sample. Thus, the 27 spectroscop- 
ically confirmed LAEs plus 62 photometrically selected 
candidates comprise our photometric sample of 89 LAEs. 
When discussing the Lya equivalent-width distribution 
of the confirmed LAEs in subsection 5.2, we use the 28 
LAEs which pass the photometric selection criteria. On 
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Fig. 8. Distribution of objects in the R — z' vs i' — TV 73816 
plane. The dots indicate objects with TV73816 < 26.0. If an 
objects is undetected (i.e., fainter than the 2a limit) in a 
given bandpass, its magnitude in that bandpass is replaced 
with the 2a magnitude. Objects on the R — z' = 1.16 line are 
those undetected in both R and z' . The red circles represent 
the spectroscopically confirmed LAEs (TV = 34); open circles 
show that they are undetected in R. The blue triangles and 
green crosses indicate nearby objects (TV = 19) and unclear 
objects (TV = 10), respectively. The yellow stars correspond 
to the 175 Galactic stars given in Gunn &; Stryker's (1983) 
spectrophotometric atlas. The black solid curves show the 
tracks of model LAEs over a redshift range of 5.5 < z < 5.9 
with four different intrinsic EW values of Lya in the rest 
frame: Wj, = 0, 20, 50, 150, and 300 A from left to right; 
filled squares correspond to z = 5.7. The magenta solid lines 
indicate the selection boundary (in the case of rrtR < jtir (2er)). 

the other hand, the contamination and completeness of 
the photometric sample are estimated using the 56 ob- 
jects lying in the high-quality region. 

The criterion of i' - TVB816 > 1.5 picks up LAEs with 
rest- frame equivalent widths of W' l > 20 A(or > 10 A). 
We adopt this criterion in order (i) not to select objects 
with very weak or no Lya emission, and (ii) to avoid a high 
contamination by foreground objects whose i' — TV 73816 
colors happen to exceed the criterion due to photometric 
errors. This lower limit to EWs in our selection is close to 
those set by Ajiki et al. (2003), W ohs > 180 A (Wo > 27 A), 
and by Rhoads & Malhotra (2001), W ohs > 75 A (W > 11 
A). However, our lower limit to EWs is slightly higher 
than that used by Hu et al. (2004), I - NB816 > 0.7 
- 1.0 depending on the R — z' color of the object. For 
LAEs at z = 5.7, Hu et al.'s criterion corresponds to 
i' — NB816 ^1.1 - 1.4. However, we find in our data 
that changing the criterion between i' — NB816 > 1.0 and 
> 1.5 does not significantly change the total number of 
candidates selected, if the selection is limited to objects 
brighter than TVB816 = 25.5; for objects fainter than this, 
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the contamination by foreground objects due to photomet- 
ric errors is expected to be high at i' — NB816 < 1.5 (see 
Fig. 7), and thus we cannot make a reliable comparison 
of the number of LAEs. For NB816 < 25.5, the condition 
i' - NB816 > 1.0 leaves 66 objects while i' - NB816 > 1.5 
selects 53. Considering that the i' — NB816 > 1.0 candi- 
dates will include a non-negligible fraction of foreground 
objects even at iVB816 < 25, the difference in the number 
of true LAEs between the two i' — NB816 limits will be 
small. 

Wc estimate the contamination and completeness of this 
photometric sample on the basis of the spectroscopic ob- 
servations. The completeness of object detection, which is 
another completeness to be considered in calculating the 
LF, is discussed in the next section. As shown in Tab. 2, 
56 out of the 63 objects with spectroscopic observation are 
from the region used to construct the photometric sample 
(725 arcmin 2 ). Among them, 36 satisfy the above crite- 
ria for LAE selection; 27 are LAEs, one nearby, and eight 
unclear objects. If all of the unclear objects are LAEs in 
reality, the contamination rate is calculated to be 1/36. 
On the contrary, if all of the unclear objects are nearby 
objects, then the contamination is (l + 8)/36. The aver- 
age of these two extremes, 14%, is regarded as an estimate 
of contamination. The completeness is estimated as fol- 
lows. Among the 56 spectroscopic objects lying in the 
region used to construct the photometric sample, 32 are 
confirmed LAEs and eight are unclear. Then, 27 out of the 
32 LAEs and all eight unclear objects satisfy the photo- 
metric selection criteria for LAEs. Thus, the completeness 
is calculated as (27 -)- 8)/ (32 + 8) if all of the unclear ob- 
jects are actually LAEs. It is reduced to 27/32 if all of the 
unclear objects are nearby objects. The simple average of 
these two values is 86%. 

These contamination and completeness estimates will 
not be so accurate, since the spectroscopic targets have 
not been chosen very uniformly in the two-color plane. 
When deriving the LFs of LAEs in the next section, we 
do not correct them for either completeness or contamina- 
tion, because the completeness and contamination correc- 
tions using the values obtained above almost cancel; the 
combined factor is (1 - 0.14)/0.86 = 1.00. 

5. Results and Discussion 

In this section, we derive the Lya LF and the far UV 
LF of z — 5.7 LAEs from the photometric sample, and 
examine the distribution of equivalent widths using the 
spectroscopic sample. The sky distribution of LAEs in the 
photometric and spectroscopic samples is also discussed. 

5.1. Lya Luminosity Function 

5.1.1. Calculation of the Luminosity Function 

Previous studies derived the Lya LF by simply divid- 
ing the observed number counts of LAE candidates in 
a given narrow band by the effective survey volume de- 
fined as the FWHM of the bandpass times the area of 
the survey. Although this procedure is accurate when the 
narrow-band filter has an ideal, boxcar shape, the shapes 



of actual filters used in LAE surveys are rather close to 
a triangle. In such cases, the narrow-band magnitude of 
LAEs of a fixed Lya luminosity varies largely as a func- 
tion of redshift. The selection function of LAEs in terms of 
equivalent width also changes with redshift; the minimum 
EW value corresponding to a given (fixed) narrow-band 
excess, such as i' — NB816, becomes larger when the red- 
shift of the object goes away from the central redshift of 
the bandpass. 

In order to avoid such uncertainties, we perform the 
following Monte Carlo simulations to find the best-fit 
Schcchtcr parameters for the Lya LF of z = 5.7 LAEs. 
First, we generate LAEs according to a given set of the 
Schechter parameters (a, 4>* , L*) uniformly in comoving 
space over the redshift range 5.6 < z < 5.8. Then, we 'ob- 
serve' these LAEs in the NB816, R, i' , and z' bands, add 
to their flux densities photon noise corresponding to the 
actual observation, and select LAEs by the same criteria 
as for selecting the actual LAEs. Thus, wc have a mock 
catalog of LAEs for the given set of (a,4>* , L*). Finally, 
we compare the NB816-band number counts derived from 
the mock catalog with the observed counts, and compute 
the x 2 value. By performing this set of simulations over a 
wide range of the Schcchtcr parameters, we find the best- 
fit parameters. 

Three points should be noted in these simulations. 
First, in the actual simulations we search for the best-fit 
</>* and L* values for only three values of a, —1,-1.5,-2, 
since the magnitude range of the observed number counts 
is not wide enough to place strong constraints on all three 
parameters simultaneously. These three a values are cho- 
sen following Malhotra & Rhoads (2004). The difference 
in x among the three a values is found to be insignificant, 
although a = — 1 gives the lowest x 2 value. We adopt the 
a = — 1.5 results as the fiducial set of the best-fit Schcchtcr 
parameters. 

Second, in order to compute flux densities in the indi- 
vidual bandpasses from a given Lya luminosity, we need 
to give an equivalent width. We assume a Gaussian distri- 
bution, /(W^dW^ozcxpi-W}) 2 /2a% v )dW i , for the prob- 
ability distribution of EWs for z = 5.7 LAEs, and generate 
objects by randomly sampling EWs from the Wq > part 
of this distribution. The value of aw can be inferred from 
the observed i' — NB816 distribution and EW distribu- 
tion. Figure 9 compares the predicted i' — NB8\Q distri- 
bution with the observation for four aw values: aw = 50 
A, 100 A, 200 A, and 400 A. For each prediction, a = —1.5 
and the best-fit L* are used. It is found that aw ~ 100 A 
and 200 A are preferred among the four values. We then 
calculate the equivalent width distribution for aw = 100 
A and 200 A, to find that the distribution for aw = 200 A 
is roughly consistent with the observed distribution of the 
spectroscopic sample, while that for aw = 100 A is much 
skewed toward lower EWs and thus clearly inconsistent 
with the observation. Based on these results, we decide 
to use aw = 200 A to find the best-fit Schechter parame- 
ters. We have confirmed that using aw = 100 A does not 
significantly change the best-fit parameters. 
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Fig. 9. Histogram of i' — ^^3816 for the photometric sam- 
ple. The four curves indicate the predictions from simulations 
assuming four values for aw'- a w = 50,100,200,400 A. 

Third, we correct the observed number counts for de- 
tection completeness, since all objects are detected in the 
simulations irrespective of their magnitudes. We esti- 
mate the detection completeness as a function of apparent 
NB816 magnitude, by distributing pscudo objects on the 
NB816 image after adding photon noise, detecting them 
by SExtractor, and computing the detection fraction. We 
assume Gaussian profiles for pscudo objects whose size 
distribution (measured by SExtractor) matches that of 
the observed LAEs. The completeness thus computed is 
> 0.8 for NB816 < 25.0, and is ~ 0.75 at NB816 = 26.0. 
The correction is thus modest in the whole magnitude 
range of the photometric sample. 

We show in Figure 10 the observed number counts by 
circles and the predicted counts by three curves corre- 
sponding to the three a values. The dotted, solid, and 
dashed lines correspond to a = —1,-1.5, and —2, respec- 
tively. All three models reproduce the observed counts 
well, although shallower a gives a slightly better fit. This 
suggests that the Lya LF of z = 5.7 LAEs is approximated 
well by the Schechter function. The best-fit parameters 



4* [Mpc~ 
< 10 42 , 



1) = (5.2 



+ 1.4 



6.3 



+3.0 



x i r>42 i O+0.4 
1.1 x 1U i L - z -0.3 



are (L* [erg s 1 ] 

10 3 ), (7.9t^2 ~ ±v J u - u -2.o 
10 43 , 1.6±J' 4 x 10- 4 ) for a = -1, 
tively. The L* values are 'observed' values. We assume 
that half of the intrinsic luminosities arc erased due to 



10- 4 ), (1.6±°;° x 
-1.5, —2.0, respec- 



(0 

a) 
u 



OA 
(0 

£ 

LO 



10 



0.1 



- 








1 


1 : 








































































































/// 










h 












// 


























5 


























a = 


-2 _ 










a = 


-1-5 










a = 

1 


-1 



22 



23 



24 25 
ATB816 



26 



Fig. 10. Number counts of LAEs in the photometric sam- 
ple. The filled circles represent the counts after correction for 
detection completeness, while the open circles show the raw 
counts. The dotted, solid, and dashed lines indicate, respec- 
tively, the best-fit Schechter functions for a = —1, —1.5, and 
—2 derived in Subsection 5.1. 

Hi gas; one can thus double L* to obtain intrinsic values. 

The best-fit Schechter functions obtained above for 
a = —1, —1.5, and —2 arc shown in Figure 11 by the thick 
dotted, solid, and dashed lines, respectively. To give some 
idea about the luminosity range in which our Schechter 
functions are reliable, we plot a simple estimate of the 
luminosity function for our sample by the filled circles. 
These data points are obtained by dividing the number 
counts in NB816 (after correction for detection complete- 
ness) by the effective volume. Here, we calculate the Lya 
luminosity of each object by subtracting the continuum 
emission measured from the z'-band magnitude from the 
total luminosity in the NB816 band, assuming z = 5.7. 
This is a rough estimate as described above, since NB816 
magnitude for a given set of L and EW varies with red- 
shift. 

From Fig. 11, the Schechter functions are found to be 
reliable over w3x 10 42 - 3 x 10 43 erg s _1 . In this range, the 
LFs for all three a values keep rising with decreasing L, 
suggesting that there exist many faint LAEs beyond our 
detection limit. Indeed, Santos et al. (2004) detected a 
high number density of faint (L ~ lfj 40-42 erg s -1 ) LAEs at 
z = 4.5 — 5.7 in a survey of lensed LAEs in very small areas. 
Malhotra & Rhoads (2004) found that a combination of 



No. ] Lyman a emitters at z = 5.7 




10 42 1 43 1 44 



L(Lya) [erg s '] 

Fig. 11. Lya luminosity functions of LAEs at z = 5.7. The 
dotted, solid, and dashed lines indicate, respectively, the 
best-fit Schechter functions for a = —1,-1.5, and —2 derived 
in Subsection 5.1. The filled circles show rough estimates 
obtained by dividing the number counts of the photometric 
sample by the effective volume. The filled triangles, crosses, 
open triangles, and open circles indicate the data given by 
Rhoads & Malhotra (2001), Ajiki et al. (2003), Ajiki ct al. 

(2004) , and Hu ct al. (2004). The dot-dashed line represents 
the luminosity function derived in Malhotra & Rhoads (2004) 
by fitting the Schechter function to a compiled data set. The 
pentagons represent the LF of LAEs at z ~ 4.5 by Dawson 

(2005) . 

Santos ct al.'s LF with brighter estimates is consistent 

with Schechter functions of a ~ — 1 2. 

Note that the simple estimate of the LF (filled circles) 
agrees fairly well with the best-fit Schechter functions. 
This means that the simple method gives a good approxi- 
mation to the true LF at least for our data. This does not, 
however, reduce the importance of the accurate derivation 
of the LF made above, since we cannot judge the reliabil- 
ity of the simple estimate without an accurate calculation. 
In addition, the simulations used to derive the LF is use- 
ful in predicting the redshift distribution of LAEs in a 
flux-limited sample as well as placing a constraint on the 
distribution of EWs. 
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5.1.2. Comparison with Previous Measurements 

Overplotted in Figure 11 are four previous results for 
z = 5.7 LAEs, which are taken from a compilation by 
Ajiki et al. (2004). The filled triangles represent the LF 
given in Rhoads & Malhotra (2001), which is based on a 
narrow-band sample of 18 candidates from a blank field 
of 710 arcmin 2 . The crosses indicate the data in Ajiki et 
al. (2003) from 20 candidates detected in a 720 arcmin 2 
area including the QSO SDSSp J104433.04-012522.2 at 
z = 5.74. The open triangles show the data in the same 
QSO field but from a survey with an intermediate-band 
filter (Ajiki et al. 2004). The open circles indicate Hu et 
al.'s (2004) measurements based on spectroscopically con- 
firmed 19 LAEs in a blank field of 702 arcmin 2 . The dot- 
dashed line denotes the LF given in Malhotra & Rhoads 
(2004) derived from a fit of the Schechter function to a 
compiled data set, including Rhoads & Malhotra (2001), 
Hu et al. (2004), and Ajiki et al. (2004). All but Ajiki et 
al. (2003) have not corrected for detection completeness. 

At L > 1 x 10 43 erg s _1 , our LF is in an acceptable 
agreement with those by Ajiki ct al. (2003) and Hu et 
al. (2004). Ajiki et al.'s (2004) measurements are lower 
than ours, probably because they adopted a higher thresh- 
old for EWs, Wq > 45 A, than the other authors includ- 
ing us. At L < 1 x 10 43 erg s -1 , our measurements are 
higher than all of the others by up to an order of mag- 
nitude. We infer that the main reason for this discrep- 
ancy is that our data are much deeper than the others 
and that the other measurements are not corrected for 
detection completeness. Large-scale inhomogeneity in the 
spatial distribution of LAEs may be another likely reason 
for the discrepancy. The differences in selection criteria 
for LAEs between our sample and the others are probably 
not a main reason for the discrepancy, since, for our data, 
changing the i' — NB816 threshold over 1.0 - 1.5 does not 
largely change the number of LAEs selected. 

The open pentagons in Figure 11 represent the LF of z ~ 
4.5 LAEs obtained by Dawson (2005) from 76 objects with 
spectroscopic redshifts. This LF is found to be slightly 
lower than our z — 5.7 LF, but the difference is within the 
la errors for most data points. Thus, drastic evolution in 
the LAE LF from z ~ 4.5 to z = 5.7 is ruled out. Ouchi 
et al. (2003) show that the LF at z = 4.9 derived from 
a sample of 87 LAEs agrees with that for z = 3.4 LAEs 
by Cowie & Hu (1998), while Maier et al. (2003) found a 
significant decline in the LAE LF from z = 3.4 to 2 = 5.7 
in their samples of z — 4.8 and z = 5.7. Our finding is 
consistent with Ouchi et al.'s (2003). The decrease found 
by Maier ct al. (2003) could be due to the small statistics 
of their samples (N = 5 and 11 for z = 4.8 and 5.7), or 
an incompleteness in object detection in their relatively 
shallow data. 

The largest uncertainty inherent in our LF measure- 
ments may be cosmic variance. The presence of large- 
scale structures of LAEs has been reported (Shimasaku 
et al. 2003, 2004; Ouchi ct al. 2005). Hu et al. (2004) 
also found a highly structured distribution of LAEs at 
z = 5.7 in their 26/5 x 26/5 field. Ouchi ct al. (2005) 
used a catalog of z = 5.7 ±0.05 LAEs for a continuous 1 
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deg 2 area to calculate the rms fluctuation of LAE over- 
density within a sphere of (comoving) 20 Mpc radius to 
be 40 ± 20 %. This is the estimate of the cosmic variance 
of LAEs on the largest scale reported so far. Our survey 
volume, 1.8 x 10 5 Mpc 3 , is five times larger than the vol- 
ume of a sphere of 20 Mpc radius. Thus, the uncertainty 
in the number density of LAEs in our sample is likely to 
be much smaller than 40%, unless there exists an unusual 
inhomogeneity in the LAE distribution on scales of 3> 20 
Mpc. 

5.2. Equivalent Widths 

5.2.1. Calculation of Equivalent Widths 

We measure the rest-frame EWs of the spectroscopi- 
cally confirmed 34 LAEs by two independent methods. 
The first method uses the z'-band magnitude to measure 
the continuum flux density, and calculates the EW by di- 
viding the continuum-subtracted flux in NB816 by the flux 
density of the continuum. The second calculates the EW 
and the continuum flux density simultaneously from the 
i' and NB816 magnitudes. Note that at z = 5.7 the Lya 
line enters both i' and NB816 bands. Malhotra & Rhoads 
(2002) used the first method, while Hu et al. (2004) ap- 
plied both and then adopted the first method. In cither 
method, we do not correct for the effect of Hi absorption 
to the Lya emission, while taking account of the IGM ab- 
sorption of the continuum emission shortward of the Lya 
line (Madau 1995). In what follows, W^(z') and 
denote respectively the rest-frame EWs calculated from 
the first and the second methods. 

The EW of the majority of our LAEs is highly uncer- 
tain because of their faint continuum emission; uncertain- 
ties are in general larger for objects with larger EWs. 
Furthermore, the probability distribution function of EW 
for each object, P(W)dW, is not symmetric but has a 
skewed shape with a long tail of large EWs, since the 
continuum emission enters into the denominator of the ex- 
pression of EW. For each object and each of the above two 
methods, we carefully evaluate P(W) using Monte Carlo 
simulations to calculate the central value of EW and its la- 
errors (see, e.g., Dawson et al. 2004). In the first method, 
we first assume that each flux density of NB&1Q and z 1 
has a Gaussian probability density function centered on 
the measured flux density with a standard deviation of 
the measured error. Then we run Monte Carlo simula- 
tions to obtain a large set of (NB816, z') flux densities. 
Finally, we calculate EWs for individual (NB816,z') com- 
binations of the set and derive P(W). Similar simulations 
are performed for NB816 and i' in the second method. 

The Wl(z') and W%{i') values of the 34 LAEs derived 
above are given in Table 3. Three values in each of the 
Wq(z') and Wg(i') columns are the central value (W c ), 
la lower limit (W_), and la upper limit (W+) from left 
to right; they are defined as f™ c P(W)dW = 0.5 and 
P(W)dW = J™* P(W)dW = 0.34. Figure 12 plots 
W${i') against Wq(z') (both are central values). The 
crosses correspond to the six objects with B or V detec- 
tion; they are identical to the six objects not passing the 
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Fig. 12. Vyg(z') versus Wg(i') for the 34 LAEs in the spec- 
troscopic sample. The crosses indicate the six objects not 
meeting the photometric selection criteria for LAEs. The 
dotted lines indicate the maximum allowable value by nor- 
mal star formation. 

photometric selection criteria for LAEs. For low EW ob- 
jects Wq(i') correlates with Wq(z'), but little correlation 
is seen at > 100 A. Some objects with W%(z') > 100 A have 
an extremely large W^i'). To derive Wp(i'), the contin- 
uum emission is calculated by subtracting the flux falling 
into the NB816 band from that falling into the i' band. 
Hence, the continuum emission calculated for objects with 
intrinsically faint continua (large EWs in most cases) can 
be close to zero or negative because of large photometric 
errors in i'-band magnitudes. Since a negative value is ar- 
tificial and it actually means an extremely large (positive) 
EW, we assign +oo to objects with negative W c . 

The relative robustness against photometric errors be- 
tween the two methods depends on the characteristics of 
the imaging data. We examine the reliability of Wfj(z') 
and VFq(i') measurements using the simulations made to 
derive the Lya LF in Subsection 5.1. Figure 13 plots 
Wg(z') (panel [a]) and W^(i') (panel [b]) against the input 
value for 2,000 model LAEs randomly selected from the 
simulated objects which satisfy the selection criteria. The 
best-fit Schechter parameters for a = —1.5 and aw = 200 
A are adopted. It is found that Wq(z') shows a broad 
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Fig. 13. Measured equivalent widths plotted against true 
(input) values for simulated LAEs. Panel (a). Wfj(z') versus 
0.5 X Wj, (input). Panel (b). W%(i') versus 0.5 X Wj, (input). 
For both panels 2,000 objects are plotted. They are se- 
lected randomly from the simulated catalog with the best-fit 
Schechter parameters for a = —1.5 and aw = 200 A. All of 
the objects satisfy the selection criteria for LAEs. 
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Fig. 14. Panel (a). Distributions of input and measured 
EW values for all simulated LAEs. The thin solid curve in- 
dicates the input (i.e., true) distribution. The dotted and 
dashed curves are for Wq(z') and WJ(i'), respectively. For 
Wg(i'), objects with Wg(i') > 500 A are shown at IV* = 505 
A. Panel (b). EW distribution of the 28 confirmed LAEs 
meeting the photometric selection criteria. The solid his- 
togram represents the number frequency of Wfj(z') (central 
values), while the solid curve corresponds to the sum of their 
probability distribution functions (P(W)). In both panels, 
the dotted line indicates the maximum allowable value by 
normal star formation. 



but tighter correlation with the input value than Wq(£'). 
Panel (a) of Figure 14 presents the relative number distri- 
bution of the input and measured EW values for all simu- 
lated objects satisfying the selection criteria. The Wq(z') 
distribution (dotted curve) agrees well with the input dis- 
tribution (thin solid curve), while the Wq(z') distribution 
(dashed curve) deviates largely in the sense that the num- 
ber density of Wq ~ 100 - 200 A objects is underestimated. 
Thus, we conclude that Wq(z') is a more reliable measure- 
ment for our data. In what follows, we adopt VFq(z') for 
the observed values. 

We have inferred that the six objects detected either in 
B or V are contaminated by a foreground object. If this is 
the case, their equivalent widths derived from NB816 and 
z 1 photometry will be lower than the true values because 
of overestimated continuum emission (the effect on Lya 
luminosity is little, as NB816 and z' flux densities are 
similarly contaminated). This underestimation of EWs 



is, however, expected to be small, since (i) among the six 
objects three are detected in z' with both B — z' > 1 and 
V — z' > 1 and (ii) three are undetected in z' but two of 
them have B and V magnitudes very close to the detec- 
tion limits. In addition, we do not use these six objects 
in the discussion of the EW distribution in the next sub- 
section. Thus, we have not corrected their EW values for 
foreground contamination. 

5.2.2. Equivalent Width Distribution 

Figure 15 plots Wq(z') against L(Lya) for the 28 con- 
firmed LAEs satisfying the photometric selection criteria. 
Objects shown by open circles have z' magnitudes fainter 
than the la magnitude. The horizontal dotted line in Fig. 
15 corresponds to a maximum value allowed by reasonable 
stellar populations (W = 240 A) derived by Malhotra 
& Rhoads (2002). Continuously star- forming populations 
with population-I or -II mctallicitics assuming a Salpeter 
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Fig. 15. Wgfz') plotted against Lya luminosity for the 28 
confirmed LAEs which meet the photometric selection crite- 
ria. Objects shown by open circles have z' magnitudes fainter 
than the lcr magnitude. The dotted line indicates the maxi- 
mum allowable value by normal star formation. 

mass function over the normal mass range cannot have 
EWs exceeding this value, unless they are as young as 
< 1 x 10 7 yr (Chariot & Fall 1993; Malhotra & Rhoads 
2002). The number frequency of W%{z') for the 28 LAEs 
is plotted in Figure 14 (b) by the solid histogram. The 
solid curve in Fig. 14 (b) corresponds to the sum of the 
probability distribution functions of the 28 objects. 

Fig. 15 (and Fig. 14 [b]) shows that LAEs at z = 5.7 
have very large EWs. The median Wq(z') among the 28 
objects is 117 A (or Wq(z') = 233 A), with lower and up- 
per quartilcs of 86 A and 147 A (171 A and 294 A in Wq), 
respectively. If the probability distribution function (solid 
curve in Fig. 14 (b)) is used, the median, lower quartile, 
and upper quartile are found to be Wq(z') = 213,125,386 
A, respectively. The EW distribution based on the spec- 
troscopic sample will be biased toward large values, since 
the spectroscopic sample appears to be biased toward red 
i' — NB816 colors (i.e., large EWs). As an independent 
check, we calculate the median and quartiles from the 
best-fit simulations with aw = 200 A, to find VP = 147 A 
(median), 80 A (1 ower quartile), and 238 A (upper quar- 
tile). Although smaller than the corresponding values 
from the spectroscopic sample, these values are still quite 
large. 



Fig. 16. Same as Fig. 15, but plotted against the far ultra- 
violet absolute magnitude, Mtjv. See Subsection 5.3 for the 
calculation of Mtjv. Objects fainter than the 2<r magnitude 
in z' are placed at Mtjv = —19.5 and shown by crosses. The 
dotted curve indicates the lower limit of Wq(z') detected in 
our sample (z = 5.7 is assumed for simplicity); objects be- 
low this curve are fainter than ./V.B816 = 26.0 and thus not 
included in our sample. 

We do not find a significant correlation in VTg(z') vs 
L(Lya) (Fig. 15). However, we find that objects with 
fainter far-UV continuum luminosities tend to have larger 
EWs. Figure 16 plots Wq(z') against far-UV absolute 
magnitude, A/tjv, f° r the spectroscopic sample. Sec 
Subsection 5.3 for the calculation of Mtjv- The dot- 
ted curve indicates the lower limit of Wq(z') detected 
in our sample; objects below this curve are fainter than 
NB816 = 26.0 and thus not included in our sample. It is 
found that all objects brighter than Mtjv = —20.7 have 
M^o( z ') — 50 A. On the other hand, the majority of ob- 
jects fainter than Mtjv = —20 have Wq(z') > 100 A unless 
a large number of objects exist below the dotted curve. 
If we interpret that objects with brighter Mtjv ar e more 
massive, then the trend seen in Fig. 16 implies that 
less massive objects have larger EWs. Recently, Ando 
et al. (2005) obtained 11 spectra of z ~ 5 LBGs with 
M uv 



<. 



-20.8, finding that bright (M< -21.5) LBGs have 
EWs less than Wq ~ 20 A while the faintest two have large 
EWs of Wq ~ 40 - 80 A. This finding is very close to the 
trend seen in Fig. 16. It is also interesting to note that a 
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Fig. 17. Same as Fig. 15, but plotted against the FWHM 
of the Lya line. The FWHM values have been corrected for 
spectral resolution. The dotted line indicates the maximum 
allowable value by normal star formation. 

similar trend has been found for LBGs at z ~ 3 although 
their absolute EW values are much lower than our LAEs' 
(Shapley et al. 2003). We also find a weak, negative cor- 
relation between Wq(z') and the FWHM of the Lya line 
(Fig. 17). This is qualitatively consistent with the trend 
seen in Fig. 16, if the line width is due to Dopplcr broad- 
ening. 

If one wants to reproduce the observed large median 
value by an ordinary, metal-enriched star-forming popu- 
lation with a Salpeter IMF (i.e. , x = 2.35) of an upper mass 
cut of 120M Q , one must assume very young ages (< 1 x 10 7 
yr), negligible dust extinction, and a nearly complete ab- 
sorption of ionizing photons in the system (Kudritzki et al. 
2000; Malhotra & Rhoads 2002). Such young ages mean 
that we are selectively picking up objects which have just 
started its star formation (or just entered a short star- 
formation phase if they have an episodic star-formation 
history). This leads us to an idea that there exist a much 
larger number of ancestors having not started star forma- 
tion and descendants with ages since the onset of star for- 
mation too old to be selected as LAEs. This idea appears 
to be inconsistent with the observed far UV LFs of LAEs 
and LBGs. In the next section we show that essentially 



all z ~ 6 LBGs pass our selection criteria for LAEs. This 
implies that nearly half of LBGs have Wq > 200 A and 
thus younger than 1 x 10 7 yr. It looks unlikely that such 
a high fraction of LBGs fall into a so narrow range of age, 
since the LBG selection picks up star-forming galaxies ir- 
respective of their ages as long as their apparent far UV 
flux densities are higher than a given detection limit. On 
the other hand, a small portion of our LAEs with equiv- 
alent widths lower than ~ 100 A can be accounted for by 
ordinary Salpeter populations of reasonably old ages as 
long as dust extinction is negligible. 

Other possibilities to account for the large median 
EW are top-heavy IMFs or zero-metallicity populations. 
Malhotra & Rhoads (2002) show that both an x = 0.5 
IMF of 1/20 solar metallicity and a Salpeter IMF of zero 
metallicity can produce EWs > 200 A even if they are 
10 s yr or older, more than ten times older than the ages 
constrained for ordinary Salpeter populations. A special 
configuration of the interstellar medium (ISM) may also 
boost the equivalent width. Hansen & Oh (2005) argue 
that Lya photons can preferentially escape from a dusty, 
multi-phase ISM than continuum photons if most of the 
dust is in cold neutral clouds, and thus that the EW of 
the transmitted spectrum can be larger than the EW of 
the unprocessed spectrum (see also Neufcld 1991). If this 
is the case, the observed large EWs can be reproduced by 
older (or steeper IMF slope) populations. 

AGNs can also produce extremely large EWs. It is, 
however, unlikely that AGNs dominate in our sample, 
since our LAEs have very narrow rest-frame widths of 
the Lya emission line, with a mean and rms of 1.07 A 
and 0.68 A, or 264 km s _1 and 168 km s _1 in velocity, 
after correction for spectral resolution. Recently, Wang et 
al. (2004) reported that no evidence of AGN activity (in- 
cluding low luminosity quasars) has been found for LAEs 
at z ~ 4.5 obtained from a narrow-band survey, either in 
deep X-ray data from Chandra or in optical spectra from 
Keck (see also Malhotra et al. 2003). Their report is con- 
sistent with our speculation, as we are probably selecting 
a similar population at z = 5.7 to their LAEs. 

5.2.3. Objects Exceeding W = 240 A 

Objects with W > 240 A are particularly interesting, 
since even a population of a very flat IMF of x = 0.5 (but 
metal enriched) cannot produce such large EWs unless it 
is younger than 1 x 10 7 yr (Malhotra & Rhoads 2002). 
We find in the confirmed 28 LAEs that 46 % (13 objects) 
have W Q (z') > 240 A. Similarly, the fraction calculated 
from the probability distribution functions of EW is 44 
%. We can also estimate the fraction from the simulations 
performed in subsection 5.1.1. The simulations for the 
best-fit Schechter parameters of the Lya LF (for a = — 1.5 
and aw = 200 A) show that 24% of the objects satisfy- 
ing the selection criteria have input values of W > 240 
A. Similarly, 29% of the objects satisfying the selection 
criteria have 'measured' values of W (z') > 240 A. These 
values are lower than the observed values, 46 % and 44 
%, probably reflecting the bias in i' — NB816 color in the 
spectroscopic target selection. Alternatively, the Simula- 
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tions could underestimate the true fraction, since we have 
not been able to constrain aw so strongly. 

From these results, we infer that the true fraction of 
W > 240 A is 30% - 40%. This value is much lower than 
that obtained by Malhotra & Rhoads (2002), 60%, from 
a photometric sample of z = 4.5 LAEs detected in their 
narrow-band survey. Our result may be closer to Hu et 
al.'s (2004); they found in their z — 5.7 LAE sample that 
the majority (11 out of 15 objects) have EWs less than 
240 A, although they did not clearly present the fraction. 
The reason for the difference between our fraction and 
Malhotra & Rhoads' is not clear. The evolutionary effect 
is unlikely, since it implies that the fraction of large-EW 
objects declines with redshift, which is opposite to the 
increasing fraction of strong Lya emitters in LBGs (see 
the next subsection). 

Recently, Dawson et al. (2004) obtained a more reli- 
able fraction of large-EW objects for z ~ 4.5 LAEs from 
a sample of 17 objects with spectroscopic confirmation, 
but they calculated the fraction of Wq > 240 A in our 
definition, not W > 240 A. Their calculation may be re- 
garded as a conservative lower limit, since Lya emission 
is assumed to be absorption free. Using the probability 
distribution function of EW from Monte Carlo simula- 
tions, they found with 90 % confidence that 18 - 29 % 
of LAEs in their sample exceed 240 A. A similar calcu- 
lation for our 28 spectroscopic LAEs finds that 7 - 25 % 
exceed W%(z') = 240 A with 90 % confidence. Thus, even 
without absorption correction, a certain fraction of z = 5.7 
LAEs cannot be accounted for by normal star formation; 
the fraction is similar to that found for z ~ 4.5 LAEs. It 
is worth noting that in our sample we do not see evidence 
of large-EW objects being free from absorption. Among 
the six objects with Wq > 150 A, at least five have a Lya 
profile which shows clear asymmetry with a sharp cutoff 
in the blue side, suggesting a large amount of absorption. 

5.3. Far- UV Luminosity Function 

5.3.1. Calculation of the Luminosity Function 

We derive the LF of the far UV continuum emission 
using the photometric sample of 89 LAE candidates. We 
transform the z'-band magnitude into the far UV contin- 
uum at the rest-frame 1350 A assuming that all of the 
LAEs are located at z = 5.7. For objects at z = 5.7, the 
z' band covers rest-frame wavelengths of 1250 A - 1480 
A with a weighted center of 1350 A. Then the far UV LF 
is calculated by dividing the number of LAEs in each 0.5 
magnitude bin by the effective volume corresponding to 
the FWHM of the NB816 bandpass (1.80 x 10 5 Mpc 3 ). In 
this calculation, the detection completeness in NB816 is 
taken into account. We do not convert the Lya LF ob- 
tained in Subsection 5.1 into the far UV LF, since the 
conversion is sensitive to the EW distribution, which is 
still uncertain for our purpose here. For reference, the 
Muv values of the 34 LAEs with spectroscopic redshifts 
arc presented in Table 3. 

Figure 18 plots the far UV LF obtained from our pho- 
tometric sample by filled circles. The vertical dotted line 
at Muv = —19.58 corresponds to the 2a limiting magni- 
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Fig. 18. Far ultraviolet (~ 1350 A) luminosity functions of 
LAEs at z = 5.7 and LBGs at z ~ 6. The filled circles represent 
our measurements from the photometric sample. The vertical 
dotted line at Mtjv = —19.58 corresponds to the 2a limiting 
magnitude in the z' band; objects fainter than this magnitude 
have not been used. The open circles show the measurements 
by Hu et al. (2004). The crosses indicate the measurements 
for LBGs at z ~ 6 given in Bouwens et al. (2005). 

tude in the z' band, 27.04 mag; objects fainter than this 
magnitude have not been used to derive the LF. Note that 
our sample is not z'-magnitudc limited but NB816-band 
limited (A^£?816 < 26.0). Since the scatter between z 1 and 
NB816 magnitudes is large, the detection completeness in 
terms of z' magnitude (or far UV magnitude) gradually 
drops when z' magnitude goes fainter because of undetec- 
tion of objects with relatively fainter NB816 magnitudes 
(for a given z 1 magnitude). We find, from the plot of ap- 
parent z' magnitude versus NB816 magnitude, that the 
measurement of the far UV LF is incomplete to some de- 
gree at Mtjv fainter than ~ —20.5. Therefore, the appar- 
ent flattening of the LF at Muv > —20.5 seen in Figure 
18 may be partly due to this incompleteness. The open 
circles show the far UV LF of z = 5.7 LAEs obtained by 
Hu et al. (2004). Their measurements agree with ours 
very well at Muv < —21, but their value at Muv = —20.3 
is lower than ours by a factor of four, probably reflecting 
the detection incompleteness of their sample. 



Lyman a emitters at z = 5.7 



17 



5.3.2. Lya Emission of LBGs 

We compare the far UV LF of our LAEs with that of 
z ~ 6 LBGs. Determinations of the far UV LF of z ~ 6 
LBGs have been recently accumulating (e.g., Dickinson et 
al. 2004; Yan & Windhorst 2004; Bunker et al. 2004; 
Bouwens et al. 2005; Shimasaku et al. 2005); see a com- 
pilation given in Bouwens et al. (2005). The differences 
in </>(M) among the determinations are roughly within a 
factor of three. We adopt here Bouwens et al.'s (2005) 
results (crosses in Fig. 18), since they are based on a 
very large sample of 506 LBGs from several deep surveys 
including the GOODS fields and the Hubble Ultra Deep 
Field (HUDF), and since their results are in the middle 
among the previous determinations. 

LBGs at z ~ 6 are detected as i'-dropout galaxies in 
multi-color broad-band surveys, and thus they include all 
galaxies with far UV continua brighter than a lower limit 
defined by the survey depth. In Fig. 18, the LAE LF is 
found to be comparable to the LBG LF over the whole 
magnitude range examined, except for our data point at 
the faintest magnitude. As stated above, the data point 
at the faintest magnitude probably suffers from a large 
incompleteness, and thus the true difference from the LBG 
LF will be smaller. We can thus conclude that at z ~ 6 
nearly all of far-UV selected galaxies with Muv ^5 — 20 
have a strong Lya emission with W > 20 A. In other 
words, an LAE survey like ours can pick out essentially 
all galaxies brighter than Muv ~ — 20 at z ~ 6. 

Malhotra et al. (2005) compared the far UV LF of 
z ~ 6 LBGs in the HUDF with the Lya LF of z = 5.7 
LAEs given in Malhotra & Rhoads (2002), to find similar 
cj>* values (after correction for an overdensity seen in the 
HUDF LBGs). They argue that the LBG space density 
could be consistent with that of LAEs, but they could not 
compare the space densities of these two populations in a 
common far UV absolute magnitude range. We demon- 
strate in this study that the space densities of LBGs and 
LAEs at z ~ 6 arc similar at least down to Muv ~ —20. 

Figure 18 may suggest that the majority of z ~ 6 LBGs 
have dust-free, extremely young ages and/or top-heavy 
IMFs, if we recall that our LAEs have extremely large 
EWs. A recent observation of far UV colors of i'-dropout 
galaxies (i.e., z ~ 6 LBGs) by Stanway et al. (2005) shows 
steep (i.e., blue) UV slopes, which might suggest young 
ages (~ 10 7 yr) or top-heavy IMFs, although their sample 
is limited to faint objects (Muv > —20) and thus does not 
overlap with our sample. On the other hand, Spitzer Space 
Telescope observations of bright (Muv ^ —21) i'-dropout 
galaxies (LBGs at z ~ 6) give old ages of a few x 10 8 yr and 
large stellar masses of a few x10 10 Mq (Yan et al. 2005; 
Eyles et al. 2005; see also Mobasher et al. 2005 for the 
finding of a J-dropout galaxy at z ~ 6.5 with an old age 
and an extremely large stellar mass) . Spitzer observations 
of HCM 6A (spectroscopically confirmed, Lya emitting 
galaxy at z = 6.56 behind Abell 370) give very young ages 
of an order of 10 6 yr (Chary, Stern, & Eisenhardt 2005), 
but this object is fainter than —21. 

It is also suggested that bright LBGs have small EWs; 
Bunker et al. (2003) and Stanway et al. (2004) made 



spectroscopy of six z ~ 6 LBGs with Afuv 1z —21, find- 
ing that two have Wq = 20 - 30 A and the rest have no 
detectable Lya emission (A similar trend was found for 
z ~ 5 LBGs by Ando et al. 2005). This supports our 
finding that LAEs with brighter Muv have on average 
smaller EWs (Fig. 16). These results collectively suggest 
that there is a large variety in the star formation history 
in z ~ 6 galaxies depending on their mass; massive galax- 
ies are considerably evolved while low-mass galaxies tend 
to be very young. 

Our study reveals that LBGs at z ~ 6 have much larger 
EWs than those at z ~ 3. Steidel et al. (2000) and Shapley 
et al. (2003) found that 20 - 25% of bright (M uv < -20) 
LBGs at z ~ 3 have W > 20 A (Their measurements in- 
clude Lya absorption). This value is a factor 4-5 lower 
than that obtained above for z ~ 6 LBGs, since almost all 
z ~ 6 LBGs pass our LAE selection, W ^ 20 A. Moreover, 
the fraction of LBGs with W > 100 A is - 2% at z ~ 3 
(Shapley et al. 2003), but this fraction rises to as high 
as about 80% at z ~ 6, since 24 out of the 28 spectro- 
scopic LAEs have W > 100 A (The probability distribu- 
tion functions of the 28 LAEs give 82%). Our finding of 
high fractions of largc-EW LBGs is qualitatively consis- 
tent with the recent finding by Nagao et al. (2005); they 
found that at least three out of 48 i'-dropout galaxies at 
z ~ 6 have W > 200 A. 

The extremely high fraction of large-EW LBGs at z ~ 6 
suggests strong evolution of Lya properties of far-UV se- 
lected star- forming galaxies over z ~ 6 and z ~ 3. The 
origin of this evolution cannot be specified in this study, 
but possible candidates will include lower dust extinction, 
younger stellar ages, and drastic changes in metallicity 
and the IMF slope toward lower and flatter values, re- 
spectively, with increasing redshift. The increase in the 
fraction of large EW objects with redshift seems qualita- 
tively in accord with the recent finding that the far UV 
continuum of LBGs is on the average bluer at z ~ 6 than 
at z ~ 3 (Bouwens et al. 2005; Stanway et al. 2005). 

We calculate the star formation rate for the 89 LAE 
candidates from the far-UV absolute magnitude using the 
formula given in Madau et al. (1998). This formula has 
been applied to LBGs and LAEs at various redshifts in 
previous studies. Our calculation is limited to the objects 
brighter than Muv = —19.58. The star formation rates 
span ~ 4 - 4OM yr -1 . By dividing the sum of the star 
formation rates by the survey volume, we obtain the star 
formation rate density to be 2.3 x 10 _3 M Q yr -1 Mpc~ 3 . 
This value is twice as high as the estimate given by Ajiki 
et al. (2003), reflecting the deeper limiting magnitude 
of our data, and thus is regarded as a new lower limit 
of the far-UV based star formation rate density in the 
z ~ 6 universe from LAE surveys. Our value is about one 
third that derived from the LF of z ~ 6 LBGs (i'-dropout 
galaxies) down to 0.3L* =3 (or —19.6 mag) by Bouwens et 
al. (2005), ~ 6 x lO~ 3 A/ yr" 1 Mpc" 3 . 

Note, however, that the formula used above to con- 
vert Muv into star-formation rate, given in Madau et al. 
(1998), assumes a constant star formation of solar metal- 
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licity with a Salpeter IMF of ages of 3> 10 7 yr. Thus, this 
formula cannot be applied to our LAEs if they have very 
young (< 10 7 yr) ages or top-heavy IMFs or zero metal- 
licity as suggested from the large EW values. The star 
formation rate will be underestimated by this formula for 
very young populations, while it will be overestimated for 
populations with top-heavy IMFs. Therefore, while we 
can robustly conclude that LAEs have a significant con- 
tribution to the total far-UV luminosity density in the 
z ~ 6 universe, the star formation rate density itself could 
contain a large systematic error. Furthermore, recalling 
that nearly all LBGs at z ~ 6 have EWs large enough to 
be selected as LAEs, we suspect that the star-formation 
rate density measurements for LBGs may also suffer from 
a similar uncertainty. 

5.3.3. Lya Bias for the LF of z ^ 6 LBGs 

Finally, let us discuss an interesting implication of the 
large EWs of z ~ 6 LBGs for their far UV LF. In most 
cases, LBGs at z ~ 6 were selected from their red i' — z' 
colors, and their Mtjv magnitudes were calculated from z'- 
band photometry (An important exception is the work by 
Shimasaku et al. 2005). At z = 6, however, the redshiftcd 
Lya line is located around the peak transmission of the z' 
band. In what follows we take the F850LP band of ACS 
as a typical z' band. If an LBG at z ~ 6 has Wq = 200 A, 
then its F850LP-band magnitude will be about 0.8 mag 
brighter than the value expected in the case of no Lya 
emission, since its observed-frame equivalent width will 
be as large as 1300 A (exact values depending on red- 
shift), which is nearly the same as the FWHM of the 
F850LP band. This simple calculation implies that the 
far-UV LF of z ~ 6 LBGs calculated from z'-band mag- 
nitudes will be biased toward higher luminosities. This 
bias is stronger for fainter LBGs, since fainter LBGs have 
larger EWs on average (Fig. 16). Comparing the LF at 
z ~ 6 with those at lower redshifts will require correction 
to the former for this bias. The correction not only will 
strengthen the recent finding that the number density of 
bright LBGs drops from z ~ 3 to z ~ 6 (e.g., Ouchi ct al. 
2004; Shimasaku et al. 2005; Bouwcns ct al. 2005), but 
also might find a new trend that faint LBGs also decrease 
in number toward z ~ 6. Nagao et al. (2004, 2005) have 
also argued significant brightening of z'-band magnitudes 
of i'-dropout galaxies due to large EWs. 

5.4- Spatial Distribution 

The sky distribution of our LAEs is shown in Figure 
19. The small dots indicate the 89 candidates in the pho- 
tometric sample. The open circles and crosses represent, 
respectively, 34 LAEs in the spectroscopic sample with 
z < 5.7 and z > 5.7. The thick solid lines outline the 
region used to make the photometric sample. Some spec- 
troscopic objects do not have a counterpart in the pho- 
tometric sample because they are either out of the region 
used to construct the photometric sample, or ruled out by 
the selection criteria. 

A large-scale density contrast is seen in the sky distri- 
bution of the 89 candidates. The average surface num- 
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Fig. 19. Sky distribution of the LAEs in the photometric 
and spectroscopic samples. The small dots indicate the 89 
candidates in the photometric sample. The open circles and 
crosses represent, respectively, 34 LAEs in the spectroscopic 
sample with z < 5.7 and z > 5.7. The thick solid lines outline 
the region used to make the photometric sample. 

ber density over the whole region is 0.12 ±0.01 arcmin -2 . 
When the whole region is split into two at X = 17', the 
surface number density is 0.097 ±0.016 arcmin -2 for the 
eastern half and 0.15±0.02 arcmin -2 for the western half; 
these values differ from the average about ±20%. The 
line-of-sight distribution of the spectroscopic sample also 
shows inhomogeneity. Objects with z < 5.7 dominate in 
the western half, while those with z > 5.7 dominate in 
the eastern half. These findings collectively suggest the 
existence of large-scale structure in our sample. 

Figure 20 plots the angular correlation function for the 
photometric sample (filled circles). Clustering signals are 
found not to be significant; the large-scale inhomogeneity 
found above is not strong enough to be clearly reflected 
in angular correlation. We also calculate the angular cor- 
relation for a bright subsample of NB816 < 25.5 (open 
circles), since brighter galaxies are in general known to be 
clustered more strongly at high redshift as well, but we 
do not see significant clustering. We, however, find possi- 
ble clustering of the ~ 2a level at a very small separation 
of < 10" in the angular correlation function, which might 
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Fig. 20. Angular correlation function of LAEs in the photo- 
metric sample. The filled and open circles correspond to the 
whole sample (N = 89) and a bright (JVB816 < 25.5) sample 
(N = 53), respectively. The inset shows an expanded view at 
> 1'. The data points have been shifted to the horizontal di- 
rection by —0.01 and +0.01 for the whole and bright samples, 
respectively, for clarity. 

suggest that LAEs tend to form close pairs. 

Ouchi et al. (2005) have examined the sky distribution 
of z — 5.7 LAEs in a 1 deg 2 area, about five times larger 
than our survey area, and found filamentary large-scale 
structures on scales of several tens of Mpc. The large- 
scale structure found in our sample appears not to be 
as prominent as the filamentary structures of Ouchi et 
al. sample. A similar observation has been obtained for 
z ~ 4.8 LAEs by Shimasaku et al. (2003, 2004). They 
have surveyed for LAEs in z = 4.79 ± 0.04 and z = 4.86 ± 
0.03 slices of the universe in a 25' x 45' area of the SDF, 
and found a large-scale structure in the latter slice, while 
finding no angular clustering in the former slice. 

6. Summary 

We have presented the results of a deep survey of Lya 
emitters (LAEs) at z = 5.7 in the Subaru Deep Field. A 
total of 89 LAE candidates with NB816 < 26.0 were se- 
lected from imaging data of six bandpasses, including the 
narrow band NB816, taken with Suprime-Cam on Subaru. 
These candidates have i' — NB816 > 1.5, red R — z' colors, 
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and B and V fluxes fainter than the 2a limits. Spectra 
of 39 objects satisfying the photometric selection crite- 
ria for LAEs were obtained with FOCAS on Subaru and 
DEIMOS on Keck II. The spectroscopic classification of 
objects was made using weighted skewness, which quan- 
tifies the asymmetry of the emission line. Among the 39 
objects, 28 are confirmed LAEs, one is a nearby galaxy, 
and eight are unclassified. Among the 28 confirmed LAEs, 
one is outside the region used to construct the photomet- 
ric sample of LAEs. Also obtained were spectra of another 
24 NB816-cxcess objects in the field; six are LAEs and 18 
are nearby galaxies. We have thus 34 confirmed LAEs 
in total. The spectroscopic results show that the photo- 
metric sample of 89 LAEs is highly reliable with a low 
contamination. 

Using these photometric and spectroscopic samples of 
LAEs, we have studied the LFs, the equivalent width dis- 
tribution, and the spatial distribution of LAEs at z = 5.7. 
We summarize the main results below. 

• By generating mock catalogs of LAEs, we have 
searched for the Schechter parameters of the Lya 
LF of z = 5.7 LAEs which best fit the observed 
number counts in the NB816 band. We have found 
(L* [ergs' 1 ], c/>* [Mpc" 3 ]) = (5.2+^ x 10 42 1 ^ +0A 
10" 3 ), (7.9± 3 ;° x 10 42 , 6.3± 3 ;° x lO" 4 ), (1.6j 
10 43 , 1.6±q; 4 x 10~ 4 ) for a = -1, -1.5, -2.0, respec- 
tively, over the range where our data are available, 
L ~ 3 x 10 42 - 3 x 10 43 erg s" 1 . The Lya LF of 
z = 5.7 LAEs is reproduced well by the Schechter 
function. Our LF is consistent with the previous 
measurements at L > 1 x 10 43 erg s~ 4 within the er- 
rors, but largely overshoots them at L < 1 x 10 43 erg 
s _1 , probably reflecting the detection incomplete- 
ness of the previous surveys. No significant evolu- 
tion is found in the LF between z ~ 4.5 and z = 5.7. 

• We have measured rest-frame Lya equivalent widths 
for the spectroscopically confirmed LAEs. We have 
found that the median value among the 28 LAEs sat- 
isfying the photometric selection criteria is W = 233 
A, with lower and upper quartiles of 171 A and 294 
A, respectively, where W , which we regard as the 
intrinsic value, is twice the measured value. About 
30% - 40% of LAEs at z = 5.7 are inferred to exceed 
W = 240 A, from both the spectroscopic sample 
and the simulations. LAEs with W > 200 A are 
unlikely to be accounted for by metal enriched, con- 
tinuously star-forming populations with a Salpeter 
IMF; metal-free populations or top-heavy IMFs are 
required. We have also found that LAEs with fainter 
far-UV magnitudes have larger EWs. 

• We have derived the far-UV (~ 1350 A) LF for 
LAEs at z = 5.7 down to Myv — —19.6 from the 
photometric sample. We have found that this LF 
agrees with the far-UV LF of LBGs at similar red- 
shifts over the whole magnitude range where the 
completeness of our sample is expected to be rel- 
atively high. This implies that at z ~ 6 almost all 
star- forming galaxies brighter than Muv ~ — 20 have 
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Lyev emission strong enough to be detected as LAEs 
(Wj, > 20 A). Moreover, the fraction of z ~ 6 LBGs 
with W\ > 100 A is found to be about 80%. These 
extremely high fractions are in sharp contrast to 
lower-redshift LBGs, suggesting drastic evolution in 
Lya emission properties. We have also argued that 
the far-UV LF and the star formation rates of star- 
forming galaxies at z ~ 6 measured from z'-band 
photometry could contain large systematic errors. 
• We have found large-scale structure in the spatial 
distribution of our sample both in the tangential and 
radial directions, although it is not as prominent 
as the filamentary structures seen in Ouchi et al.'s 
(2005) much larger survey for z = 5.7 LAEs. 
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ID 


a (J2000.0) 


6 (J2000.0) z I 


Sw 


NB81Q B 


V 


R i' 


z' i' - NB&16 R - z' 


(1) 


(2) 


(3) (4) (5) 


(6) 


(7) (8) 


(9) 


(10) (11) 


(12) (13) (14) 



18699 
20087 
20495 
23759 
27787 
31858 
31765 
36334 
39849 
42576 
46904 
48328 
50215 
61418 
70600 
73078 
84305 
84720 
93966 
95588 
96007 
108164 
113271 
121315 
124530 
124783 
126848 
132343 
152586 
154296 
158036 
163079 
166310 
168127 



13 24 38.94 
13 24 40.52 
13 24 11.88 
13 24 7.20 
13 23 59.71 
13 24 56.80 
13 24 15.99 
13 25 23.41 
13 23 59.17 
13 24 5.27 
13 24 16.47 
13 23 55.36 
13 24 11.87 
13 24 24.11 
13 23 54.60 
13 25 13.27 
13 24 6.39 
13 23 58.98 
13 25 5.06 
13 24 30.63 
13 24 33.09 
13 24 0.15 
13 24 23.70 
13 24 37.19 
13 25 6.48 
13 25 22.12 
13 24 45.48 
13 25 15.60 
13 24 11.89 
13 24 13.00 
13 24 0.51 
13 23 51.96 
13 24 16.13 
13 23 59.27 



+27 13 40.9 
+27 13 57.9 
+27 14 2.0 
+27 14 41.6 
+27 15 26.3 
+27 16 9.7 
+27 16 11.1 
+27 17 1.3 
+27 17 40.7 
+27 18 11.6 
+27 19 7.7 
+27 19 23.7 
+27 19 48.2 
+27 22 15.1 
+27 24 12.7 
+27 24 42.0 
+27 27 4.2 
+27 27 13.4 
+27 29 10.3 
+27 29 34.5 
+27 29 38.6 
+27 32 12.1 
+27 33 24.8 
+27 35 2.3 
+27 35 44.5 
+27 35 46.8 
+27 36 12.6 
+27 37 19.9 
+27 41 31.8 
+27 41 45.8 
+27 42 35.3 
+27 43 38.2 
+27 44 11.6 
+27 44 33.8 



5.645 
5.724 
5.697 
5.687 
5.707 
5.692 
5.691 
5.705 
5.692 
5.684 
5.665 
5.693 
5.682 
5.748 
5.654 
5.700 
5.695 
5.682 
5.721 
5.738 
5.696 
5.672 
5.710 
5.716 
5.689 
5.720 
5.718 
5.697 
5.681 
5.715 
5.676 
5.676 
5.698 
5.674 



D 
D 
D 
F 
F 
D 
F 
D 
F 
F 
F 
F 
F 
D 
F 
D 
F 
F 
F 
F 
F 
F 
F 
F 
D 
D 
F 
D 
F 
F 
F 
F 
F 
F 



5.03 
8.18 
7.03 
5.81 
6.03 
6.01 
5.90 
3.75 
4.85 
6.30 

10.30 
7.28 
6.48 
8.00 
8.77 

11.62 
8.38 

20.66 
5.36 
8.60 
6.28 
4.13 

11.82 
3.38 
3.33 
7.66 

10.49 
6.28 
6.11 
5.01 

12.68 
4.38 

10.61 
6.28 



±0.82 
±0.23 
±0.41 
±5.16 
±2.25 
±0.88 
±2.52 
±0.43 
±1.14 
±0.98 
±1.65 
±4.88 
±1.94 
±0.01 
±0.80 
±2.14 
±1.84 
±6.02 
±3.74 
±4.55 
±5.63 
±2.23 
±1.01 
±2.33 
±0.19 
±0.80 
±6.73 
±0.83 
±1.64 
±2.76 
±1.53 
±1.21 
±1.80 
±0.78 



25.80 
24.00 
24.77 
24.65 
25.23 
25.57 
24.38 
24.62 
24.21 
24.54 
24.15 
25.80 
25.30 
25.00 
24.81 
25.79 
24.28 
24.47 
25.18 
25.41 
25.20 
25.89 
23.41 
25.37 
24.62 
25.20 
25.35 
24.78 
24.18 
25.57 
24.83 
24.54 
23.93 
24.27 



28.8 
28.8 
28.8 
28.8 
28.8 
28.8 
28.8 
27.6 
28.8 
28.8 
28.8 
28.7 
28.8 
28.8 
28.8 
28.8 
28.8 
28.8 
28.8 
28.3 
28.7 
28.8 
28.8 
28.8 
27.4 
28.8 
28.8 
28.8 
28.8 
28.8 
28.8 
28.8 
26.1 
28.8 



28.1 
28.1 
28.1 
28.1 
28.1 
28.1 
28.1 
26.9 
28.1 
28.1 
28.1 
28.1 
28.1 
28.1 
28.1 
28.1 
28.1 
28.1 
28.1 
28.1 
28.1 
28.1 
28.1 
28.1 



28.20 
28.20 
28.20 
28.20 
28.20 
28.20 
28.20 
26.47 
28.20 
28.20 
28.20 
28.20 
28.20 
28.20 
28.20 
28.20 
28.20 
28.20 
28.20 
28.05 
28.20 
28.20 
27.91 
28.20 
27.5 27.17 
28.1 28.20 
28.20 
28.20 
28.20 
28.20 
28.20 
28.20 
25.83 
28.20 



28.1 
28.1 
28.1 
28.1 
28.1 
28.1 
25.9 
28.1 



27.84 
26.48 
27.52 
27.40 
27.84 
27.84 
26.96 
25.93 
27.38 
27.45 
26.18 
27.49 
27.84 
27.83 
26.70 
27.84 
26.96 
27.45 
27.84 
26.93 
27.47 
27.84 
25.80 
27.77 
27.01 
27.20 
27.71 
27.61 
27.03 
27.84 
27.84 
27.16 
25.43 
27.06 



27.04 
26.39 
27.04 
26.26 
27.04 
27.04 
26.53 
25.34 
26.82 
26.87 
25.37 
27.04 
27.04 
27.00 
27.04 
27.04 
26.49 
26.94 
27.04 
26.79 
27.04 
27.04 
24.73 
27.04 
27.04 
26.28 
27.04 
27.04 
26.63 
27.04 
27.04 
26.87 
24.84 
27.04 



1.84 
2.22 
2.33 
2.24 
2.31 
2.21 
2.13 
1.17 
2.86 
2.67 
1.69 
1.63 
2.36 
2.24 
1.72 
1.95 
2.23 
2.17 
2.58 
1.21 
1.72 
1.87 
2.10 
2.04 
1.88 
1.75 
2.19 
2.38 
2.24 
1.86 
2.48 
2.11 
1.19 
2.34 



1.16 
1.81 
1.16 
1.94 
1.16 
1.16 
1.67 
1.13 
1.37 
1.33 
2.83 
1.16 
1.16 
1.20 
1.16 
1.16 
1.71 
1.26 
1.16 
1.26 
1.16 
1.16 
3.17 
1.16 
0.13 
1.92 
1.16 
1.16 
1.57 
1.16 
1.16 
1.33 
0.99 
1.16 



(b) nearby objects 



ID 


a 1 


;j2ooo.o) 


8 (J2000.0) 


z 


I 


28247 


13 


24 


11.70 


+27 


15 


30.5 




F 


29275 


13 


25 


17.89 


+27 


15 


46.4 




D 


31925 


13 


25 


15.61 


+27 


16 


11.0 


1.182 


D 


34775 


13 


25 


8.10 


+27 


16 


48.1 




D 


38133 


13 


23 


56.28 


+27 


17 


26.1 


1.180 


F 


42561 


13 


24 


3.41 


+27 


18 


17.8 


0.630 


F 


56181 


13 


25 


25.87 


+27 


21 


12.0 


0.620 


D 


59788 


13 


25 


10.24 


+27 


21 


53.2 


0.628 


D 


68251 


13 


23 


40.67 


+27 


23 


45.6 


0.635 


D 


70071 


13 


24 


34.91 


+27 


24 


10.2 


0.623 


D 


76702 


13 


24 


5.75 


+27 


25 


37.3 


1.186 


D 


78892 


13 


24 


15.26 


+27 


25 


60.0 


0.615 


F 


96705 


13 


24 


25.59 


+27 


29 


47.2 


0.632 


F 


99588 


13 


23 


57.78 


+27 


30 


30.1 


0.636 


F 


110439 


13 


25 


24.72 


+27 


32 


44.3 


0.629 


D 


122518 


13 


25 


13.45 


+27 


35 


17.9 


1.174 


D 


136295 


13 


25 


5.54 


+27 


38 


10.0 


0.244 


D 


165225 


13 


24 


53.38 


+27 


43 


57.7 


0.637 


D 


168136 


13 


24 


3.01 


+27 


44 


34.7 




F 



Sv 



NB816 B V 



R 



-NB8W R 



2.99 ±0.30 
-6.40 + 0.12 
-8.49 + 0.20 
0.14±0.19 
-1.88 ±0.05 
0.65 ±0.01 
1.68 ±0.23 
1.34 ±0.14 
1.44 ±0.00 
0.26 ±0.00 
-5.99 ±0.00 
-0.23 ±0.01 
2.92 ±0.67 
-1.82 ±0.43 
0.06 ±0.00 
-1.03 ±0.00 
-3.23 ±0.01 
-13.95 ±0.77 
-3.69 ±0.33 



24.05 
22.79 
24.40 
23.57 
22.23 
21.83 
24.88 
25.02 
23.46 
22.35 
22.93 
23.10 
25.27 
24.21 
23.45 
25.35 
22.27 
24.67 
23.88 



27.0 
25.8 
26.4 
28.1 
24.0 
24.1 
27.5 
27.8 
25.8 
24.6 
24.6 
25.8 
28.5 
26.9 
25.9 
28.8 
24.2 
27.0 
25.5 



26.9 26.83 
25.4 24.87 
26.1 25.88 
27.4 26.60 
23.8 23.73 

23.8 23.41 
27.4 26.79 

28.0 27.05 
25.6 25.04 
24.4 24.04 
24.4 24.39 

25.6 24.95 

27.9 27.89 
26.8 26.41 
25.8 25.36 

28.1 28.20 
24.0 23.69 

26.7 26.10 
25.4 25.30 



26.21 
24.35 
25.62 
25.55 
23.51 
23.14 
26.23 
26.68 
24.73 
23.77 
24.17 
24.42 
27.04 
25.99 
24.99 
27.84 
23.71 
25.81 
25.09 



26.82 


2.00 


0.01 


23.52 


1.02 


1.35 


25.17 


1.02 


0.71 


24.52 


1.18 


2.08 


23.21 


1.03 


0.52 


23.22 


1.02 


0.19 


26.54 


1.03 


0.25 


27.04 


1.50 


0.01 


24.87 


1.10 


0.17 


24.08 


1.24 


-0.04 


23.96 


1.04 


0.43 


24.93 


1.08 


0.02 


27.04 


1.69 


0.85 


26.34 


1.42 


0.08 


25.46 


1.41 


-0.10 


27.04 


1.89 


1.16 


23.89 


1.30 


-0.20 


25.83 


1.04 


0.26 


24.79 


1.03 


0.51 



22 



Shimasaku et al. 



(c) unclear objects 



ID 


a 


J2000.0) 


5 (J2000.0) z 


I 


Sw 


NB816 


B 


V 


R 


i> 


z' %' 


-NB816 R-z 1 


17721 


13 


25 


9.51 


+27 


13 


29.4 ... 


D 


2.93 + 0.53 


24.67 


28.8 


28 


1 


28.09 


26.81 


25.45 


1 


73 


2.63 


34503 


13 


25 


7.24 


+27 


16 


38.7 ... 


D 


0.92 + 1.30 


25.50 


28.8 


28 


1 


28.20 


27.84 


26.84 


1 


94 


1.36 


41967 


13 


23 


58.35 


+27 


18 


1.4 ... 


F 


0.01 + 0.28 


25.73 


28.8 


28 


1 


28.20 


27.84 


27.04 


2 


03 


1.16 


74901 


13 


23 


58.83 


+27 


25 


5.1 ... 


F 


-2.86 + 2.14 


25.49 


28.8 


28 


1 


28.20 


27.49 


27.04 


1 


88 


1.16 


89624 


13 


24 


24.80 


+27 


28 


11.9 ... 


F 


2.04 + 0.97 


25.70 


28.8 


28 


1 


28.20 


27.84 


27.04 


1 


84 


1.16 


91179 


13 


24 


17.84 


+27 


28 


33.4 ... 


F 


2.20 + 0.92 


25.05 


28.8 


28 


1 


28.20 


27.84 


27.04 


2 


32 


1.16 


97631 


13 


24 


0.34 


+27 


30 


1.0 ... 


F 


2.01 + 1.40 


24.98 


28.8 


28 


1 


28.20 


27.16 


26.73 


1 


87 


1.47 


98040 


13 


24 


32.88 


+27 


30 


8.8 ... 


F 


-6.53 + 2.32 


24.22 


28.8 


28 


1 


28.20 


26.44 


25.45 


1 


83 


2.75 


105644 


13 


23 


49.49 


+27 


31 


42.7 ... 


F 


-0.13 + 0.80 


24.74 


28.8 


28 


1 


28.20 


27.05 


26.04 


1 


56 


2.16 


111685 


13 


24 


8.75 


+27 


32 


58.8 ... 


F 


-0.60 + 0.25 


25.70 


28.8 


28 


1 


28.20 


27.84 


27.04 


2 


08 


1.16 



Notes. 

(1): ID in the NB816-dctected catalog. 

(4) : Rcdshift. For LAEs 1215.67 A is adopted as the Lya wavelength. 

(5) : Instrument used. D = DEIMOS, F = FOCAS. The nearby object 76702 was observed with FOCAS as well. Wc 
have adopted the DEIMOS spectrum with a better spectral resolution. 

(6) : Weighted skewness. 

(7) - (12): Magnitudes in the six bandpasses. NB816 magnitude is the MAG_AUTO magnitude, and the others are 
2"-aperture magnitudes. Magnitudes fainter than the 2a magnitudes have been replaced with the 2a magnitudes. 
(13), (14): Colors measured on a 2" aperture. 



Table 2. Breakdown of the spectroscopic sample in terms of the selection criteria for our photometric sample 



class 


selection criteria 


total 




satisfy 


not satisfy 




LAE 


28 (27) 


6(5) 


34 (32) 


nearby 


1(1) 


18 (15) 


19 (16) 


unclear 


10 (8) 


0(0) 


10 (8) 


total 


39 (36) 


24 (20) 


63 (56) 



Notes. 

Numbers in parentheses are the number of objects lying in the region used to construct the photometric sample of 
LAEs (See Fig. 19 and section 4). 



No. ] Lyman a emitters at z = 5.7 23 



Table 3. Lya and continuum properties of confirmed LAEs 



ID z FWHM L(Lya) sp L(Lya) W%{z') Wffi) M uv SFR 
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81 


/"'A 

60 


1 O 1 
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1 o 
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84720 


5.682 


1 A A 

10.9 
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■-T.-T 
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48 
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69 


' > o 

32 
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8.4 
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C TO 1 
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3.3 
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AO 

92 


A TCT 
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oc 
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40 
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24 
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16 
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A O 
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49 


on 

29 
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24 


11 


61 
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5.8 
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1.3 
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3.3 ±1.3 


1 A T 
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68 
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1 AA 
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21 


00 






113271 


5.710 


9.5 


7.0 


25.4 ± 0.8 


39 


36 


42 


57 


49 


67 


—22.17 


40.1 


121315 


5.716 


0.0 


3.0 


4.9 ±0.8 


148 


71 


505 


52 


26 


276 






124530 


5.689 


2.2 


2.7 


9.4±1.1 


131 


78 


307 


31 


20 


61 


-19.65 


3.9 


124783 


5.720 


12.0 


4.0 


5.2 ±0.8 


35 


26 


48 


35 


22 


65 


-20.58 


9.2 


126848 


5.718 


18.8 


4.9 


5.0 ±0.8 


91 


53 


221 


73 


33 


592 






132343 


5.697 


3.6 


3.4 


8.1 ±1.0 


144 


77 


387 


98 


43 


'OC' 






152586 


5.681 


7.6 


2.6 


13.9 ±1.3 


90 


63 


146 


82 


45 


197 


-20.59 


9.3 


154296 


5.715 


4.9 


0.8 


4.0 ±0.9 


131 


60 


431 
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35 


oc 
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5.676 


11.2 


1.9 
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92 


54 
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211 


52 


00 
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5.4 
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5.676 


4.8 


1.8 


10.6 ±1.4 


96 


64 


178 


66 


35 


193 


-20.24 


6.8 


166310 


5.698 


9.3 


3.4 


12.6 ±0.9 


21 


19 


24 


8 


6 


9 


-22.08 


36.7 


168127 


5.674 


5.5 


3.9 


15.4 ±1.4 


261 


143 


656 


146 


69 


1423 







Notes. 

(1) : ID in the NB816-detected catalog. Shown in bold font are objects not satisfying the photometric selection criteria. 

(2) : rcdshift. 

(3) : FWHM of the Lya line in the observed frame (A). 

(4) : Lya luminosity (xlO 42 erg s _1 ) calculated from the spectrum. 

(5) : Lya luminosity (xlO 42 erg s _1 ) calculated from the NB816- and z'-band photometry combined with the spectro- 
scopic redshift. 

(6) : Wq(z') (A). The first, second, and third columns are the central value (W c ), lc lower limit (W-), and la 
upper limit (W+). These values are defined as f™ c P(W)dW = 0.5 and P(W)dW = P(W)dW = 0.34, where 
P(W)dW is the probability distribution of EW derived from Monte Carlo simulations. 

(7) : W^{i') (A). The first, second, and third columns are the central value, la lower limit, and la upper limit. The 
definition of these values is the same as for Wq(z'). 

(8) : Far-UV continuum absolute magnitude (AB). 

(9) : Star formation rate calculated from Muv (M Q yr _1 ). 



